Harmonic Gravitational Wave Spectra of Cosmic String Loops in the Galaxy 
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A new candidate source of gravitational radiation is described: the nearly-perfect harmonic series 
from individual loops of cosmic string. It is argued that theories with light cosmic strings give 
rise to a population of numerous long-lived stable loops, many of which cluster gravitationally in 
galaxy halos along with the dark matter. Each cosmic string loop produces a spectrum of discrete 
frequencies in a nearly perfect harmonic series, a fundamental mode and its integer multiples. The 
gravitational wave signal from cosmic string loops in our Galactic halo is analyzed numerically and 
it is found that the for light strings, the nearest loops typically produce strong signals which stand 
out above confusion noise from Galactic binaries. The total population of cosmic string loops in the 
Milky Way also produces a broad signal that acts as a confusion noise. Both signals are enhanced 
by the clustering of loops gravitationally bound to the Galaxy, which significantly decreases the 
average distance from the solar system to the nearest loop. Numerical estimates indicate that for 
dimensionless string tension G/i/c 2 < 10 -11 , many loops are likely to be found in the Galactic halo. 



Lighter strings, down to G/i/c 



10' 



are detectable by the Laser Interferometer Space Antenna 



(LISA). For these light strings, the fundamental and low-order harmonics of typical loops often lie 
in the band where LISA is sensitive, 0.1 to 100 mHz. The harmonic nature of the cosmic string loop 
modes leaves a distinct spectral signature different from any other known source of gravitational 
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I. INTRODUCTION 

Topological defects produced during cosmic phase 
transitions are a standard component of many field and 
brane theories and cosmologies. Often they take the form 
of macroscopically extended classical one dimensional ob- 
jects, with microscopic radius, called cosmic strings. The 
cosmological evolution of strings forms a population of 
quasi-stable loops that lose energy mainly by gravita- 



tional waves [fl 
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36j . This paper describes a new way 
in which these gravitational waves may be observed: a 
regime in which radiation from individual loops can ap- 
pear as perfect harmonic series in an observed frequency 
spectrum. 

String loops oscillate, radiate gravitational waves and 
shrink until they completely decay. The center of mass 
speed of a loop when it forms is of the order of unity. 
After it stops interconnecting with the rest of the string 
network, a loop's velocity decays inversely proportional 
to the scale factor. At late times, the primordial velocity 
is negligible and the loop population clusters in almost 
the same way as the dominant cold, collisionless dark 
matter (3?j]- In the case of light strings, for which loops 
are small and numerous, a galaxy halo can contain a very 



large number of loops [38| • They are concentrated to high 
density, so the mean distance to the nearest loop is much 
smaller than the cosmic mean. 

A cosmic string loop produces a spectrum of dis- 
crete frequencies [l], [39j which may be detectable if it 
is close enough. The spectrum of any loop is given by a 
sum over a nearly perfect harmonic series of frequencies 
/„ = 2nc/L, where L is the length of the string loop. 
This distinctive property is unlike any other astrophysi- 
cal source of gravitational waves and if observed, would 
provide convincing evidence of the existence of cosmic 
strings, as well as detailed information about their astro- 
physical behavior. 

The power from each discrete mode n of a loop is given 
by, 



E n = P n Gfi 2 C, 



(1) 
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where fj, is the mass (energy) density of the string. We 
define J — ^2 P n - Numerical simulations indicate 7 is ap- 
proximately 50 to 100. We use a value of 50 for this study. 
The power in each mode depends upon the particular os- 
cillation pattern of a string loop, but the general solution 
of a sum of harmonic modes does not depend upon any 
particular model. For the illustrative estimates in this 
paper we assume a very simple model, again motivated 
by numerical estimates: the mean power of an ensem- 
ble scales as (P n ) oc n~ 4 / 3 '. Loops clustered around the 
Milky Way halo, most of which are too distant to detect 
individually, taken together create an unresolved back- 
ground of gravitational waves, akin to the white dwarf 
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binaries [4CJ, |41| . We also estimate this confusion back- 
ground. 

The new effects are most important for cosmic string 
loops with lower string tensions Gp,/c 2 < 10 -12 [42L l43l| 
and for stable loops formed at a significant fraction, 
a = 0.1, of the horizon [42|, [43|, [44|. In this situation, we 
find that cosmic strings are an important new source at 
the frequencies 0.1 Hz to 10 mHz band where Laser In- 
terferometer Space Antenna (LISA) [45|,[46[ will be most 
sensitive. The extension into this regime is all the more 
relevant given the current limits set on cosmic string ten- 
sions via SDSS, WMAP, and millisecond pulars [47ll48j|. 
Since the results depend on the size distribution of stable 
loops [H, [44[ , a number of possibilities are considered in 
the parameter studies here. 



II. DISTRIBUTION OF LOOPS 

The average number density of loops between size L 
and L + dL is n(L)alL, computed numerically. The loops 
are assumed to be clustered around galaxies in the same 
way as the dark matter. The number density in the 
galaxy is matched to the dark matter halo, give n b y the 
NFW density distribution p NF w(r) 0, mM, M, H3 . 
The number density is then a function of the length of 
loops and distance from the center of the galaxy n(I/, r). 



1. Density in the Galaxy 

The density of galactic string loops is matched to the 
dark matter distribution in the Milky Way using the 
NFW density pnfw{t)'- 



PNFw{r) 



x(l +x) 2 ' 
r 



x = — , 



(5) 
(6) 



where p s and r s are determined by observation. Repre- 
sentative values are given by [50| using the favored model, 
e.g. r s =21.5 kpc. From above we see that for r << r s , 
p <x r^ 1 and for r >> r s , p oc r~ 3 . 

From the previous example using strings of tension 
Gp = 10~ 12 , N = 10 4 , r t /r s = 10, and r s = 21.5 kpc, we 
find at the earth's distance from the galactic center, r = 8 



kpc, a value of n « 9.4 x 10 kpc , orn 



-1/3 



2.2 kpc, 



an "average" distance to the nearest loop. The nearst 
loop is likely closer as the distribution is not constant 
over this range of position. 

The number density of loops in the Milky Way is 
strongly influenced by the size and mass density of the 
cosmic string loops. Tables IH- MII give results for the 
number of loops in the Milky Way for various param- 
eterizations of the string loops. The results clearly show 
the large predominance of loops for very light and large 
strings. 



A. Loop Density 

We start by figuring out the size of the loops at the 
fundamental frequency: 



L = 2c/ h 



(2) 



Using the one-scale model [5J, [55| , loops form with size 
L w acH^ 1 ^) and start to decay at a rate L = —-fGp. 
Loops created at a time t c will be of the size L(t c , t) at the 
time t. At the present time t we find the fundamental 
frequency given by: 



h 



2c 



acH-^Q-jGpito-tJ/c 



(3) 



The number density at t a , for time of creation t c is given 
by: 

Using the time t c we solve numerically to find the number 
density. 

To illustrate with typical numerical values: the average 
distance between galaxies is 5 Mpc, the number of loops 
within the galaxy for Gp = 10~ 12 and a fundamental 
frequency of 1 mHz is approximately 10 . 



III. POWER FROM EACH LOOP 



A loop radiates power in each mode n modeled by, 

E n {r') (x ^, (7) 

where E n is power radiated per mode from a loop. The 
total power is 



E = ^Gp 2 = p nGp 2 



(8) 



where 7 is found to be between 50 and 100, and the P n are 
power coefficients of the individual modes. Note that the 
power radiated is independent of the loop size. To first 
order we ignore directionality of the emitted gravitational 
radiation. 

Each loop has a spectrum given by a set of power co- 
efficients, P n , and the average sum for the population 
(%2Pn) ~ 7 is fixed by the statistical properties of the 
strings. This is an average over a wide array of loops: 
individual loops vary from this. For most loops the most 
power is in the fundamental mode and the lowest modes, 
while the higher modes have significantly reduced power. 

For a given mode the quality factor is given by, 



Ann 

PnGp 



(9) 
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For light strings with G[i = 10 12 we find, 

Qn « 10 12 ^. 



Pn 



(10) 



For the fundamental mode P n is of order unity at most, 
so Q is large, and increases with higher n. Thus the 
harmonic series is almost perfect and the lines extremely 
narrow. This justifies the use of delta functions in other 
applications. 

Given the extremely large value of Q n the signal from 
each mode in the loop is practically constant for the du- 
ration of observation with respect to decay. Potential fre- 
quency shifts due to path differences from gravitational 
lcnsing, gravitational doppler effect, Newtonian acceler- 
ation in the Galaxy, and relative motion have not been 
analyzed [30| . 



IV. STRAIN PRODUCED BY 
GRAVITATIONAL RADIATION 

The relative weakness of gravitational waves from cos- 
mic strings allows the use of linearized gravity. Assuming 
plane wave solutions, we find for the energy flux: 



F, 



- — f z h z 



(11) 



where h n is the root-mean-square strain, and using c = 1. 
Inserting c we find : 



h Jp±Ge±. 

7T C Z rf n 



h n = 3.095 x HT 12 -/^ [GM C = 1)] 

mm- 



(12) 



(13) 



For G/j, = 10 12 , f=l mHz, and r=l kpc, we find hi 
10 -21 , within the detection limits of LISA. 



V. SPECTRUM OF A SINGLE LOOP 

Using previous results for the strain produced by a 
gravitational wave, we can add the harmonic time de- 
pendence to find the strain at each mode for the result: 



h n (t) 



cVK Gjj.1 
7r/„ c 2 r 



-i2%f n t 



(14) 



Again, the P n are the "power coefficients" estimated to 
scale in the mean as n -4 / 3 . Recall that the sum of co- 
efficients is given by ^ P n — 7. In practice most of the 
power is in the fundamental mode, so the the power from 
lower modes drops off, allowing truncation of the sum at 
a reasonable value of n. As an example, we set teh am- 
plitudes equal to the sample mean: 



Pn = Pm 



-4/3 



(15) 



and Pi = 18. At n — 20 we find ^ P n « 45, giving a 
reasonable estimate. For a single loop the total strain 
will be a sum of all the modes, 



h{t) = J2h n {t). 



(16) 



In practice a suitable cutoff is used to facilitate numerical 
computations. 



Calculation of Single Loop Spectrum using 
Discrete Fourier Transform 



Since the signal detected by LISA will be as a discrete 
set of points, it is useful to take a discrete Fourier trans- 
form of a sample signal. In general LISA has a sam- 
pling rate of f sa mp = 10Hz, which leads to large num- 
bers of data points for integration times of one to three 
years. This large number of data points also leads to a 
very clear signal due to the harmonic nature of the grav- 
itaional waves produced by the cosmic strings. 

Given the signal h(j) where j is an integer and the j th 
data point corresponding to h(t) = h(jAT), we define 
the discrete-time Fourier transform as [57j |. 



h{k) 



1 



E 

3=1 



. (j-i)(fc-i) 



(17) 



N is the total number of data points and k is the integer 
value in Fourier frequency space. Note the introduction 
of the square root on the integration time T, which is 
standard in the sensitivity curves for LISA [561 ]. 

The gravitational wave signal from the cosmic strings 
has the form, 



h(j) = £ 



K.f, 



— ^-sin(2^ W iV), 



(18) 



where oj n is not the physical frequency, but the discrete 
time frequency, which is different than the physical fre- 
quency given by /„ = 2nc/L, and P n is the magnitude 
of the n th power mode. Using P n = Pin -4 / 3 this is then 
inserted into Eq. [T71 

Investigations of the parameter space of the loops in- 
dicates the heaviest individually detectable loops are of 
tension G\i w 10~ 10 for a = 0.1, due to their reduced 
numbers in the Milky Way, see Table [U Shown in Fig.[T] 
is the largest likely signal from an individual loop at the 
frequency 10 Hz. This signal stands above the noise 
from binaries. 

For the large loops, a — 0.1, it is found that the lightest 
individually detectable loops have tension G/i > 10~ 16 . 

The average size a of the loops as a fraction of the 
horizon when they formed affects the number of loops 
currently radiating. For large a the number of loops 
remaining is large, thus at a given frequency the distance 
to the nearest loop is shorter. This is indicated in Fig. [5] 
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FIG. 1: Plot of the cosmic string loop spectrum for G/i = 
10 -11 , a— 0.1, at a distance of r = 7 kpc. Observation time 
is T = 1 year, the sampling rate is 0.1 Hz, Pi = 18, and the 
fundamental frequemcy is /i = 10 -4 Hz. The first 10 modes 
are shown. This signal stands above the confusion noise from 
galactic binaries which is shown on the graph. 
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FIG. 2: Plot of two cosmic string loop spectra with G/i = 
10 -12 : the "thin" spectrum a=0.1 at a distance r = 2.2 kpc 
from the solar system, and the "thick" spectrum a = 10 _o 
with r=8.9 kpc. 



Note the amplitude varies as which shows only a 
small difference on our log plots. The total number of 
loops is significantly larger for large a. 

The increase in the radiated power of the heavier string 
loops results in greatest variation in signals potentially 
detectable. Fig. [3] shows the effect for string of tension 
G[i = 10~ 16 . The heavier string loops are, on average, 
farther away but their larger output makes up for in- 
creased distance. 

The signal to noise for a loops is given by the standard 
formula for a periodic source [591 ] . 



SWiT 



\ho\ 



T 



Sn(fY 



(19) 



where T is the integration time and S n (f) is the noise 
spectral density of the detector. For LISA the curve for 
S n (f) is given in 

As an example we find for G/i = 10~ 12 , a = 0.1, and 



log(5,,(f ) 1/2 Hz~ 1 ' 2 
-14r 




logff) Hz 



FIG. 3: Plot of the Fourier spectrum of cosmic string loops 
with tension Gfi = 10 -16 and a distance of r = 0.066 kpc. 
The fundamental frequemcy is fi = 1 mHz and the first 10 
modes are shown. Note the heavier loops (Gfj, = 10 -12 ) have 
a much larger signal, in spite of their greater distance from 
the solar system. The confusion noise from galactic binaries 
is also shown. At this fundamental frequency, the loop is not 
detectable, but shifted to higher frequencies it is. 



h = 10~ 3 that the SNR w 10 5 



VI. TOTAL GALACTIC SIGNAL 

Another key ingredient is the gravitational wave signal 
from loops within the galactic halo. These loops are not 
resolvable on an individual basis, but do contribute a 
significant gravitational wave signal. The isotropic and 
stationary background from both current and evaporated 
loops has been calculated previously. To this background 
we add the signal from loops within the dark matter halo 
of the milky way, which is not isotropic from the solar 
system. 

We sum over distances in bins of Ar' from the earth 
out to the edge of the galaxy. The total flux F net of grav- 
itational wave energy at a given frequency from strings 
of frequency fk = 2k /L in the k th mode, within volume 
dV a distance r' from the Earth is given by, 



= jFr(r')n( X ',L) dV. 



(20) 



Here we use k to label the modes to avoid confusion with 
the number density n(L). 

From the total flux at a given frequency we find the 
strain, 



fn 



7TC° 



(21) 



which should be interpreted as the rms strain at a fre- 
quency /„. 

Results of numerical calculations are plotted in 
Figs. H 03 and El 
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FIG. 4: Plot of the cosmic string loop strain spectrum for 
large loops a = 0.1 in the galaxy. The top curve is of string 
tension G/i = 10~ 12 and the bottom curve is Gfi = 1CP 20 , 
in increments of 10 2 . Also included are the LISA sensitivity 
curve with an integration time of 1 year, and the galactic 
white dwarf noise. For each loop only the fundamental mode 
is included. 



log(h) 




FIG. 5: Plot of the cosmic string loop strain spectrum for 
small loops a = 10 in the galaxy. The top curve is of string 
tension Gfi = 10 -12 and the bottom curve is G/i = 10 -20 , 
in increments of 10 2 . Also included are the LISA sensitivity 
curve with an integration time of 1 year, and the galactic 
white dwarf noise. For each loop only the fundamental mode 
is included. 



VII. CONCLUSION 

These results suggest a new way to observe light cosmic 
strings. They indicate that individual cosmic string loops 
are detectable due to the local concentration of loops 
with Galactic dark matter. They display a unique and 



distinctive harmonic spectrum which requires no special 
template fitting to detect, only a Fourier transform of the 
signal. The large quality factor (of order 1/G/x) ensures 
the frequency spectrum is almost static over the several 
years of observation, further increasing detectability of 
the signal. 

This study has taken a somewhat simplified approach 
to the distribution of the loops, simply matching them 
to the dark matter halo of the Milky Way. The model 
of the loops is also oversimplified but adequate for the 

log(h) 




FIG. 6: Plot of the cosmic string loop strain spectrum for 
small and large loops in the galaxy. It is assumed 90% of 
the loops shatter into small loops of a — 10~ 5 . The top 
curve is of string tension G/x = 10 -12 and the bottom curve 
is Gfj, = 10 -20 , in increments of 10 2 . Also included are the 
LISA sensitivity curve with an integration time of 1 year, 
and the galactic white dwarf noise. For each loop only the 
fundamental mode is included. 



purpose of estimating detectability. 

For the total galactic signal, we estimate that Galac- 
tic string backgrounds are detectable by LISA down to 
Gfj, = 1CU 19 . This is significantly more sensitive than 
previous results from the exragalactic stochastic back- 
ground. Even though the broader cosmological impact 
of the loops is minimal, they still may provide a definite 
gravitational wave signal of new physics. 
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TABLE I: Number of Loops in the Milky Way for the given 
fundamental frequencies and a— 0.1 





/(Hz) IO" 4 


10~ 3 
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10" 11 
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io- 12 


lxlO 4 


lxlO 4 


lxlO 4 


lxlO 4 


io- 13 


3 xlO 5 


3 xlO 5 


3 xlO 5 


3 xlO 5 


1Q -14 


9 xlO 6 
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lxlO 7 


10~ 15 


lxlO 8 


3 xlO 8 


3 xlO 8 


3 xlO 8 


io- 16 


3 xlO 8 


4 xlO 9 


9 xlO 9 


lxlO 10 



TABLE II: Number of Loops in the Milky Way for the given 
fundamental frequencies, varying a and G/i = 10~ 12 
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/(Hz) IO" 4 
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10~ 5 
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2 xlO 2 


2 xlO 2 
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lxlO 2 



TABLE III: Number of Loops in the Milky Way for the given 
fundamental frequencies, varying a and G/i = 10~ 16 
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/(Hz) IO" 4 


10~ 3 
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3 xlO 8 
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lxlO 10 


10~ 2 


lxlO 8 


lxlO 9 


3 xlO 9 


4xl0 9 


IO -3 


3 xlO 7 


4 xlO 8 


9 xlO 8 


lxlO 9 


IO" 4 


1 xlO 7 


lxlO 8 


3 x10 s 


3xl0 8 


10~ 5 


3 xlO 6 


4xl0 7 


9xl0 7 


1 xlO 8 


10~ 6 


1 xlO 6 


lxlO 7 


3xl0 7 


3xl0 7 



